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Abstract 

The microscopic spectral correlations of the Dirac operator in Yang-Mills theories coupled to 
fermions in (2 + 1) dimensions can be related to three universality classes of Random Matrix Theory. 
In the microscopic limit the Orthogonal Ensemble (OE) corresponds to a theory with 2 colors and 
fermions in the fundamental representation and the Symplectic Ensemble (SE) corresponds to an 
arbitrary number of colors and fermions in the adjoint representation. Using a new method of 
Widom, we derive an expression for the two scalar kernels which through quaternion determinants 
give all spectral correlation functions in the Gaussian Orthogonal Ensemble (GOE) and in the the 
Gaussian Symplectic Ensemble (GSE) with all fermion masses equal to zero. The result for the 
GOE is valid for an arbitrary number of fermions while for the GSE we have results for an even 
number of fermions. 



NBI-HE-00-19 



hep-th/0004081 



1 Introduction 



Random Matrix Theory has successfully been used to extract information about the spectral cor- 
relations of the Euclidean Dirac operator^ = 7^(5^ + iA^) eigenvalues in the low energy limit of 
Yang-Mills theories such as QCD. In (3 + 1) dimensions and in the low energy limit the effective Yang- 
Mills partition function coincides with the partition function defined by the chiral Random Matrix 
Theory (^RMT) when the so-called microscopic limit is taken |l], ^, |3], |j. The low-energy partition 
function describes the fermion mass dependence in the static limit and in a finite volume of space-time 
V [|| and is determined alone by global symmetries. The finite volume implies that we restrict to 
the case where only the low-lying excitations (the Goldstone modes) contribute to the field theory 
partition function, while the kinetic terms of the Lagrangian are neglected. In the microscopic limit 
we look at the Dirac operator spectra on the scale A = 0(V~ l ), which corresponds to a magnification 
of the spectra in the vicinity of A = on the scale V . In case of spontaneous chiral symmetry break- 
ing, reflected in a condensate different from zero, S / 0, the scale V~ l equals the average eigenvalue 
spacing. This follows from the Banks-Casher relation §1 £ = irlimV^ 1 p(ty, where p(0) is the spectral 
density of the Dirac operator evaluated in origin and where first the thermodynamic and subsequently 
the chiral limit is taken. Thus we see that spontaneous chiral symmetry breaking is intimately related 
to the spectrum of the Dirac operator in origin. 

In Random Matrix Theory no dynamical information is incorporated, only the global symmetries of 
the physical system are being used and therefore one studies the "universal" spectral correlations of 
the eigenvalues of the considered operator. Based on the work of Leutwyler and Smilga @ it has 
been conjectured that the spectrum of the Dirac operator in QCD and similar theories is universal in 
the microscopic limit jl], (8|. This conjecture is supported by the fact that the sum rules derived by 
Leutwyler and Smilga, which involve inverse powers of the Dirac operator eigenvalues, can be derived 
from RMT @, g. Thus, in this limit the microscopic spectral correlation functions (of which the 
microscopic spectral density is the most simple) of the Dirac operator can be derived from a much 
simpler Random Matrix Theory in which only the symmetries of the Dirac operator are the inputs. 
Depending on the representation of the gauge group SU(N C ) and the number of colors -/V c , the Dirac 
operator belongs to one of three universality classes, which in xRMT are represented by the orthogonal 
(xOE), the unitary (xUE) and the symplectic (xSE) ensemble ||. For each theory the symmetries of 
the Dirac operator specifies one of these universality classes. In addition the chiral structure of the 
Dirac operator in all three theories is incorporated in xRMT, requiring a specific block structure of 
the matrices in the three ensembles xUE, xOE and xSE. See ||, |9| for a more detailed discussion. 



Each of the three Yang-Mills theories in (3 + 1) dimensions has an analogue in (2 + 1) dimensions 
and the effective partition function of each theory has been showed to coincide with the partition 
function in each of the three universality classes defined by non-x Random Matrix Theory (non-^RMT) 
|4], K], 11]. In an odd number of space-time dimensions chiral symmetry does not exist. But in (2 + 1) 



dimensions with an even number of flavors Nf it has been suggested that the spontaneous breakdown 
of flavor symmetry occurs, and this is the analogue of spontaneous chiral symmetry breaking in (3 + 1) 
dimensions [f|, 12|. Thus, the argumentation of the entire picture in (3+1) dimensions has a parallel in 

11] . Here we learn again that an order parameter £ of the flavor symmetry 



§0 



(2 + 1) dimensions 

breaking in (2 + 1) dimensions is related to the spectral density of the Dirac operator, evaluated at 
zero, through a generalization of the Banks-Casher relation. The lack of chiral symmetry in (2 + 1) 
dimensions is in the three ensembles non-xOE, non-xUE and non-^SE reflected in the lack of the 
chiral block structure of the matrices. 
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Thus in both (3 + 1) and (2 + 1) dimensions we have three types of Yang-Mills theories, defined by 
a choice of fermion colors N c and a representation of the gauge group SU(N C ), and the symmetries 
of the Dirac operator in each theory implies a specific structure of the Dirac matrix. Each theory is 
represented by a random matrix ensemble: in (3+1) dimensions we have the three chiral (x) ensembles 
and in (2 + 1) dimensions the three non-chiral (non-x) ensembles. 



GAUGE GROUP 


REPS. 


DIRAC MATRIX 


ENSEMBLE 





SU(2) 


Fundamental 


Real 


(non-)xOE 


1 


SU(iV c > 3) 


Fundamental 


Complex 


(non-)xUE 


2 


SU(A C ) 


Adjoint 


Quaternion real 


(non-)xSE 


4 



Table 1: The table illustrates the classification, each ensemble is labelled by the Dyson index = 1,2, 4. 

In all three Yang-Mills theories in (3 + 1) dimensions, as well as in (2 + 1) dimensions we can have an 
arbitrary number of flavors Nf of course. However, the coincidence of the ensembles (3 = 1 and (3 = 4 
in (2 + 1) dimensions with the corresponding effective field theory is only for even Nf. The relations 



to non-xGOE and non-^GSE of the two field theories was very recently derived by Magnea |lC], 11]. 
In a theory with fermions in the fundamental representation of SU (2) the following flavor symmetry 
breaking pattern causes the creation of the condensate £ IT 



Sp(2N f ) -► Sp{N f ) x Sp(N f ). (1) 



In [11] the flavor symmetry breaking pattern 

0(N f )^0(N f /2)xO(N f /2), (2) 

is shown to be the one in a theory with N c arbitrary and the adjoint representation of SU(N C ). In 
this paper we want to derive the massless microscopic spectral density of the Dirac operator in these 
two field theories from the Random Matrix Theory of non-^OE and non-xSE. 

The universality |i~3| , 14] of the ensembles (3 = 1, (3 = 2 and (3 = 4 allows the choice of a Gaussian 
distribution (G) in these three ensembles, which is an advantage in view of the calculation of the 
spectral correlation functions. With the help of orthonormal polynomials all the microscopic spectral 
correlation functions have been derived in xGOE, xGUE and xGSE with massive fermions, see review 



in [15]. In non-xRMT, however, only the microscopic spectral correlation functions in non-xUE with 



an arbitrary number of massive fermions have been derived [|], [TJ| 17]. The two remaining universality 
classes are the non-^OE and non-^SE. In this paper we derive the kernels sffi (x, y) for non-xGOE and 
non-xGSE with massless fermions. They determine all massless spectral correlation functions in these 
two ensembles. Specifically we derive the massless microscopic spectral density in the two ensembles. 
A direct verification of our results is possible through the generation of matrices distributed according 
to the probability distribution in the ensembles non-xGOE and non-xGSE. We also compare with the 



spectral sum rules recently derived in [10, 11]. 



The traditional method to derive spectral correlation functions in general orthogonal and symplectic 
ensembles with the use of polynomials is known as Dyson's quaternion matrix method [18|. In this 
method the kernels S^\x,y), (3 = 1,4, are represented by special sums involving sfteiu-orthonormal 
polynomials and the spectral correlation functions are determined by the quaternion determinant of 
a quaternion matrix given by S^\x,y), (3 = 1,4. Widom and Tracy 24] have modified Dysons 
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quaternion matrix method, in the sense that the polynomials in the relevant kernels now can be 
chosen arbitrary. In |H| the relevant kernels for the ensembles (3 = 1 and (3 = 4 are given by 
orthonormal polynomials. The main content of section 3 in this work is to provide a simple recipe 
for how to derive the two kernels in an ensemble defined by a general weight function, only by the 
use of orthonormal polynomials. We avoid the actual proof and focus on the construction of a 
helpful machinery to derive the needed kernels. In section 4 we use the recipe on the non-xGaussian 
ensembles (3 = 1,4, defined in section 2, and thus derive an expression for the two kernels S^\x,y). 
Although the equivalences between the field theory partition function and the the partition functions 
in the two cases (3 = 1,4, only are valid for even Nf, and thus the result only has interest for even 



Nf in Yang- Mills theory in (2 + 1) dimensions |10| , 11 1, the method gives us the spectral correlation 
functions also for odd Nf in the ensemble (3 = 1. In section 5 we present results for the microscopic 
spectral densities in the two ensembles and compare them with spectral sum rules and with Monte 
Carlo simulations done directly on random matrices. 



2 The non-xGaussian ensembles 

The non-x random matrix model is defined by the partition function [Q, 10, 11 1 

N f 

zf f {M) = DT P#(D = / DT Y[ det(*T + ro/)e -*F ^ W), (g) 

f=i 

where T is from an ensemble of hermitian N x N matrices and the integration is taken over the 
Haar measure DT. The Dyson index (3 has the value (3 = 1 for the orthogonal ensemble (non-xOE), 
(3 = 2 for the unitary ensemble (non-xUE) and (3 = 4 for the symplectic ensemble (non-xSE), which 
corresponds to real, complex and quaternion real matrix elements respectively. The matrix model is 
for a generic potential V(T), but the basic assumption of universality justify the use of a Gaussian 
distribution, consistent with no additional input but the symmetries of the system. The condensate 
is E 0, the parameter Nf is restricted to integers and the diagonalized mass matrix j\A is having 
Nf masses in the diagonal. The name "non-x" is attached due to the inclusion of the determinant in 
(||), which makes the integrand not always positive and which makes these ensembles rather different 
from the usual orthogonal and symplectic ensembles near zero. 

Putting all fermion masses to zero and deleting the determinant term (Nf = 0) makes the ensemble 
equivalent to the well known (classical) Gaussian ensemble and the replacement of the determinant 
with its absolute value (or restriction to even Nf) gives the generalized Gaussian ensemble f20|| . The 
matrices in the chiral (x) ensembles are rectangular in general and have a specific block structure, 
the former a result of the incorporation of the analogue of topological charge in the Random Matrix 
Theory and the latter corresponds to the choice of a representation of the Dirac matrix in a chiral 
basis. The lack of topological charge and chiral transformations in an odd number of space-time 
dimensions is in the non-x RMT model (|3|) reflected in the quadratic matrices with no additional 
constraints, but the symmetries in each of the ensembles (3 = 1,2,4. This is exactly what separates 
the non-x ensemble from the x ensemble. 

By decomposition of the matrices T one can easily transform to integration over the N eigenvalues 
{X k } of T. Choosing the Gaussian distribution then gives the partition function 



r oo N N f/ 2 roc N 

Z n]{M) = X\d\ k \{ (Xl + m})*"?* \A({X k })f = / J[ d\ k J^(X k )\A({X k })f, (4) 
J ~°°k=i f=i J -°°k=i 



3 



for even Nf, and for odd Nf we have an extra mass m which by the choice m = gives 

r oc N (Nf-l)/2 N 

zf f {M)= l[dX k X k J] (Xl + mj)e-— x l\A({X k })f= ]JdX k uf(X k )\A({X k })f, (5) 

J -°°k=i f=i J -°° •- • 

were we have neglected all irrelevant overall factors. Here the function 

N 



k=l 



A({A 4 })^n^- A ^ 

i<j 



(6) 



is the Vandermonde determinant and we have defined the weight function of the non-x Gaussian 
ensemble 



„ N/3T. .2 

mj)e * A fc 



n (Ag + 

(iV>-l)/2 

A fc J] (Al + mJJe-^ 



for A^j even, 
for A^j odd, 



(7) 



on the support I =] — oo, oof. In the next section we outline the general method to derive the m-point 
spectral correlation function in a general ensemble defined by a weight function u. Note that, in the 
massless case, i.e all mf = 0, there is no separation between the two cases of even and odd Nf. Our 
goal is to derive an important function which determines the m-point spectral correlation functions 
in the two cases (3 = 1,4 of the massless non-^Gaussian ensembles. 



3 The m-point correlation function 

The spectrum of the ensemble defined by the general weight function u> on the support I is described 
by the m-point correlation function 

1 AM f N 
R$(xi,...,x m ) = {/3) _' / dx m+1 . . . dx N Y[iv(xj) \A({ Xi })f. (8) 
Z Nf v )■ z j=1 

This function gives the probability density that m of the eigenvalues, irrespective of their ordering, 
are located in infinitesimal neighborhoods of an, ... , x m . 

Now we will review how to derive the function R$ with the help of orthonormal polynomials in all 
three ensembles (3 = 1,2,4. Define the functions 

<Pi{x) = Pl {x)ij{x) l l 2 , i = 0,1,2,..., (9) 

where {pi(x)} is the sequence of polynomials orthonormal with respect to the weight function oo{x) 
on I. Thus the sequence {<fi(x)} consists of orthonormal functions on X. Define the Hilbert space TL 
spanned by the functions (pi(x), i = 0, 1, 2, . . . , (N — 1). The projection operator K onto the space TC 
has the kernel 

K N (x,y) ee |> fc (*V fc (|/) = ^(^W VN-i(x)) ( J J ) ( JZ^y) ) ' (10) 
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where the last equality follows from the important Christoffel-Darboux (see fl21| ) formula and cin = 
fc/V-i/&/v> k?{ denoting the highest coefficient in p^{x). In the unitary ensemble (J3 = 2) the use of 



orthonormality and a simple rewriting of |A({Aj})| gives [18] 



R$(xi,...,x m ) = det[K N (xi,Xj)] 1<ij<m . (11 



<i,j<m' 

In the orthogonal and symplectic ensembles we have an analogous result for R { £\ = 1,4. Here, 
however, the corresponding kernel KjS is not a scalar, but a quaternion kernel and Rjn is represented 



by the quaternion determinant of quaternion matrix [Kn(xi, yj)]i<i,j<m [ IS | . Now, representing the 



quaternion kernels by their 2x2 matrix representations, whose entries all are given by one specific 
scalar kernel , completes the analogous picture for the ensembles = 1,4 : The matrix kernels 

KW( xv )-i( S N&y) S^D(x,y)\ 
2 V IS y N '{x,y) S y N J {y,x) J 

and 

K W (xv) -( 4\*>V) SWD(x,y)\ 

\ IS y N '{x,y) -e(x-y) Sy(y,x) J 

determine the m-point correlation function 

R^( Xl ,...,x m ) = QdetiK^ixuXj)]^^, = 1,4. (14) 

Thus, the m-point correlation function is represented by the quaternion determinant of a 2m x 2m 
matrix. Here e(x) = x/(2\x\) (the kernel of the operator e), the scalar kernel S^\x,y) is given by 
certain sums of products involving the functions (f n , and I and D are integration and differentiation 
operators, respectively. The operator has kernel S^\x,y), S^D{x,y) is the kernel of the 
operator S^D and IS^\x,y) is the kernel of ISjfp. We see that the entire matrix kernel K^\x,y) 
is given by the scalar kernel sffl(x,y), = 1,4. In the classical way derived by Dyson p^ j, and 
later generalized by Mehta and Mahoux |f23|, the polynomials in the functions <p n , and therefore in 



S^\x,y), are chosen sfcew-orthonormal with respect to u, which through an important theorem by 
Dyson gives the result (|i~4|). In p4|], it has been shown that any choice of a family of polynomials 
leads to the same matrix kernel K^\x,y), = 1,4. The choice of skew-orthonormal polynomials 

leads apparently to the most simple Kj^\x,y). But the tedious work to derive skew-orthonormal 
polynomials and the lack of an relevant formula for this analogue of the Christoffel-Darboux formula 
(which is preferable when the scaling limit N — > oo is taken) asks for a representation of S^\x,y) 
in which the polynomials have well known properties and the sums involving them are more easy to 



deal with. This is precisely the result of [19]. For a weight function u, with respect to which there 
exist orthonormal polynomials, and for which the function u'/u is a rational function, the two kernels 
S%\x,y), 0=1,4, are given by @ : 



2n 

S$(x,y) = K 2N (x,y)- £ [A^C^^e^y), (15) 

i>n,j=l 

In 

S$(x,y)=K N (x,y)- £ [AC (I — BAC)~ 1 ]jitp i (x)£t/jj(y). (16) 

i<n,j=l 
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Both scalar kernels S^\x,y), (3 = 1,4, equal a modified version of the scalar kernel ( |i~0| ) from the 
corresponding unitary ensemble, plus a linear combination of n functions. The number n equals 
the sum of orders of the poles of uj' /uj in the extended complex plane, where every endpoint of the 
support Z, where uj' juj is analytic, is included as a simple pole. Thus, for the Gaussian ensemble 
(cj(x) = e~ x , 2 =] — oo,oo[) we have n = 1 because of the simple pole in oo and for the Legendre 
ensemble (uj(x) = 1, X = [—1, 1]) n = 2 because of the two endpoints. For the weight function uj the 
polynomials in (|l6| ) (the ensemble (3 = 1) are orthonormal with respect to uj 2 , meaning we must make 
the shift uj 1 ! 2 — > uj in (||). In the function Ql5D, the involved polynomials are orthonormal with respect 
to uj, as in the ensemble with (3 = 2, but here the shift ./V — > 2N must be taken. The relevant Hilbert 
space for each ensemble is thus 

(3 = 1 : H = span{(p ,(pi,...,(pN-i}, where <pi(x) = pi(x)uj(x), 

(3 = 4: : H = span{^o, Vl) • • • > P2N-l}, where ipi(x) = pi(x)uj(x) 1/2 , (17) 

and Pi{x) is a polynomial of degree i. The functions ipi{x), « = 1,2,..., 2n, are determined from the 
two functions </?jv(x), <pn-i{x) (N — ► 2N for (3 = 4 and uj 1 / 2 — > uj for (3 = 1), and the poles of uj' /uj. 
Once we have determined the 2n linear independent functions, that is precisely ipi £ TC, i = 1, 2, . . . , n, 
and ipi S TL- 1 , i = (n + 1), . . . , 2n, the matrices A, B and C (and Aq, Coo and Co, which just are block 
parts of A and C ) can be derived from these. In the next section we describe the procedure in great 
detail. 

The results ( |i~5| ) and ( |l~6| ) have their origin in identities for each operator (3 = 1,4, which are 

valid for a general uj, telling that equals the projection operator K onto the Hilbert space H, plus 
a correction. This correction is in each ensemble = 1,4 given by the operators K, D, and e, and in 
case the operator has finite rank, i.e independent of N, the kernel of this correction constitute 

finitely many terms. When uj' juj is a rational function this is precisely the case, and the choice of 
orthonormal polynomials with respect to uj and uj 2 for (3 = 4 and (3 = 1 respectively leads to (|l5| ) and 

(ED- 

In the proceeding section we construct a recipe for how to derive the n corrections on the right hand 
sides of (15) and ([llf). 



3.1 Deriving the kernels S^\xi, Xj), (3 = 1,4, with the help of orthonormal polyno- 
mials 



For a semi-classical weight function uj, i.e uj' juj is a rational function, the following procedure leads 
to the derivs 
on the form 



to the derivation of the two associated kernels S^\xi,Xj), (3 = 1,4. We restrict to weight functions 



fulfilling 



u;{x) = e~ v{x) , (18) 

lim w (x) = 0, (19) 
x—>ol 

where V{x) is continuously differentiable in the interior of I. (For certain weight functions the con- 
dition ( |j~9|) is unnecessary [p6f| .) The n correction terms in each of the two ensembles are determined 
by the poles of uj' juj and the two functions 

(3 = 1 : <pi(x)=jH{x)u{x), i = (N-l),N (20) 

(3 = 4 : <pi(x) = Pi (x)u{x)3, i = (2N-l),2N, (21) 
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where the polynomials pi are chosen orthonormal with respect to uj 2 (x) = e -217 ^ and uj{x) = e~ v ^ x \ 
respectively. The following steps give the right hand side of ([l^) and (|l^), where the specific functions 
given above, belonging to each ensemble, are used. In the ensemble (3 = 1 we assume that iV is even. 
With the purpose to work with both ensembles at the same time we make the shift 2N — ► N for 
(3 = 4. Therefore, in the end of the procedure we must remember to double up the matrix dimension 
N — > 2N for /3 = 4. Here is the recipe summarizing the general method [19] for obtaining the two 
kernels S$(x,y), (3 = 1,4 : 

(1) At first the main contribution to S^\x,y), f3 = 1, 4, is determined. For each ensemble = 1,4, 
the kernel Kjy(x,y) is constructed from the associated two functions, ( ^0[ ) or (21), through (10). 

(2) Using the definition 

c( « flfor^ = 4 
x - y 12 for (3 = 1, 

calculate the four functions 

Ajv(x) = a N dy <p N (x)tp N -i(y)U(x,y), (22) 



jBjv(sc) = aiv y^iiy ip 2 N (y)U(x,y), (23) 
= iiv(x) - —V'(x), C(x) = -J^-Bn^x), (24) 

l ClN-i 

for each (3 = 1, 4. (The c^ =1 ) = 2 comes from the factor 2 in the exponent of the weight function 
uj 2 .) The kernel of the operator is [D, K](x,y) and it is given by [] 



g(a)-C(y) A(x)-A(y) 



[D,K](x,y) = a N ( m (x) ^(x)) I A(g ^ (y) g^^Z^ ) ( J. (26) 

Here the central matrix is vital for the derivation of the specific matrix A (not to be confused 
with the function A(x)) later on. 

(3) Next, the number n = rioo + ^fZi n x i is calculated. The orders {n Xj } of the N p poles {xj} of the 
function <J /uj are added. A pole at infinity with order is included as well as the endpoints of 
X in which uj is analytic, the latter are counted as simple poles (n Xj = 1). (Notice that (uj 2 )' /uj 2 , 
relevant for (3 = 1, has the same poles with same orders as uj' /uj.) 

(4) By writing out the power series of the rational functions A(x), Bn(x) and C(x) in ((2^) one can 
easily verify that the kernel [D, K](x,y) can be represented by a linear combination of the 2n 
functions 

x k ip N -i(x), x k ip N (x), (0<k< noo ), (27) 



1 This is easily derived with the help of : 

<Pn(. x ) \-( A ( x ) B N (x) \ f <p N (x) , 

A-iW J - { ~C(x) -A{x) ) \ VN -!{X) I 

and the representation (|l0|) of Kn(x,u). 
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(x-Xi) k Viv-i(ar), (x-Xi) k Viv(z), (0 < A; < n*,). (28) 

In addition, the space, Q, spanned by these 2n functions has a subspace of dimension n contained 
in H and another subspace of dimension n contained in H . Therefore, we now determine In 
linearly independent functions ipi(x) G TC, i = 1, 2, . . . , n, and ipi(x) G Ti.- 1 , i = (n + 1), . . . , 2n. 
It is easily shown that all the functions in Q are orthogonal to all functions contained in the 
space 7i/S, where 

S = span{(/?7v_ fc , ipi}, < < rioo, < k < (n - n^), (29) 

is a subspace of Ti and has dimension n. This implies that the n functions ipi G Ti 1 - are given 
by the demand ipi G simplifying the derivation of the functions. 
We have 

In 

[D,K](x,y) = J2 Airfiixtyiy), (30) 

where the In x In matrix A now must be derived through (^) by the specific choice of the 
functions ipi(x), i = 1,2,..., 2n. The matrix A is symmetric (because K is symmetric and D is 
anti-symmetric) and satisfies 

Aij = if i,j<n or i, j > n, (31) 
which reduces the work. There are not other restrictions on the choice of the functions ipi, but 



the linear independence, and the matrix A is uniquely determined through (26) by the given 
choice. 

(5) The 2n x 2n matrix B defined by the inner product 



Bij = (eipi,ipj) = ( J^dye{x - y)ipi(y) , ipj) = J i J j dx d V £ ( x ~ y)V'i(y)^jOr), (32) 
must now be calculated. The fact that e{x) = —e{—x) implies 

Bij = —Bji, Ba = 0, (33) 

meaning that we only have to determine n(2n — 1) matrix elements. Note that the elements are 
numbers depending on the parameters ./V and j3. In general the result for S^\x, y) in the scaling 
limit N — > oo has the highest interest. If an asymptotic relation for the orthonormal polynomials 
(and thus for ipi) is well known, it is therefore relevant to examine whether or not it is legal 
to put in the asymptotic relations when this limit is taken in (32). It is certainly preferable to 
interchange the integration and the limit N — > oo in (32), the justification of which is given by 
Lebesgue's Majorant Theorem. The general procedure is valid for all (even) N, though, and in 
the next steps we work with the general matrix B, with no reference to the actual calculation 
of the matrix elements. 



(6) The matrices on the right hand sides of ( |l5| ) and (16) are given by the 2n x 2n matrices J and 
C 

Jij = Sij - (<5f/)|i,j> n , 0<i,j<2n, (34) 

and 

C = J + BA. (35) 
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For instance for n = 1 we have 

'-(;!)■ <»> 

The matrix Aq is defined by the 2n x n matrix produced by deleting the last n columns in A 
and the nx2n matrix Cq is defined by the matrix produced by deleting the last n rows in C. In 
addition we have the n x n matrix Coo given by deleting the last n rows and the last n columns 
in C. We need only the n last rows of the matrix 

AqCqqCq, (37) 

and the n first columns of 

AC{I - BACy 1 . (38) 

(7) The functions eipi(x) = Jjs(x — y)i^i{y) are ingredients in the n correction terms in (|i"5| ) and in 
(|l6|), as well as in the matrix elements B^j = (eipi,ipj). If allowed, it is highly preferable to put 
in the asymptotic relation for ipi(x) in the scaling limit N — > oo. We must solve the integrals 
explicitly if this is not allowed. 



(8) The results of points (1), (4), (6) and (7) are collected to construct the right hand side of ( |i5|) 
and (|i~6|). For j3 = 4 we must remember the replacement N — > 2N everywhere, giving the right 
hand side of (|T5|). 

4 The kernels S^j -(x^Xj) 



In this section we use the recipe of last section to derive the two kernels ( |15| ) and (16) in the the 
massless non-% Gaussian ensemble, defined by the choice irif = in the weight function (0). 
For all rrif = we have the weight function 

J£ ) f (x) = x N fe- cx \ = 1,4, (39) 

on the support I =] — oo,oo[, where c = "^F ■ The spectral density p{x) is the 1-point correlation 
function and from (|8|) we have P] 

= (-!)"%(*). (40) 

We will always take even making the spectral density an even function, also when Nf is odd. For 
even Nf the polynomials [2C] 

rrNf, \ , ^m/T^f.-l 1+N f T^o^-t 2^ 7^f ^(^2~ + m +\) /„ n N 
H 2ra ( X = "! ) (hm f ) ^ 4 L m 2 (cx ), ^' = -12 2/ 41 

m! 



fl&iW = (-ir(^)^c^ x L?V), ^ = r ^ + (42) 



are orthonormal with respect to (|39| ) onI=] — oo, oof. Here the functions L"(x) are the generalized 
Laguerre polynomials (which are orthonormal with respect to u>(x) = x a e~ x , a > —1, on X =]0, oo[), 
with the co efficient of x n given by 

k n = {-l) n /n\. (43) 

We call the polynomials (f4l|) and ( |42] ) generalized Hermite polynomials. Notice that these polynomials 
satisfy 

H^(-x) = H^(x), H% +l (-x) = -H% +1 (x), (44) 
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and 

<'(*) = i<^ 2 (x), (45) 

-N f 

the latter property is easily derived from a basic relation between Laguerre polynomials. Letting k i 1 
be the coefficient of the highest power in H i (x) (and remembering that N is even) implies 

/ i \ N ~* , -AT. 1 JV-2 S+^f 

fciV-l 2 ^(fcpcVc-^ ,JV .X 1 

a ^v = -r — = — -jv - N . i ~ = (it 2 c 2 , 46 

kN (-l)T fc„ (4')" 3 cf c ^ 2 

where fej^ is the coefficient (^3[). The generalized Hermite polynomials are relevant for all three cases 



(3 = 1,2,4 and even Nf. For /3 = 2, 4 and even iVy we must use the polynomials (41) and ( |42| ) which are 
orthonormal with respect to (|39"|), but for odd Nf we cannot use the orthonormal polynomial method 
described earlier, due to the fact that there does not exist orthonormal polynomials with respect to an 
odd weight function defined on an even interval [21]. For (3 = 1, however, we must use the orthonormal 
polynomials with respect to u^~ l \x) 2 = x 2N f e~ 2cx2 , implying that we can choose Nf both even and 



odd in this case. The relevant polynomials are given (41) and (|42j ) with the replacements Nf — > 2Nf 
and c — > 2c, which we will carry out at the end of derivations. Thus, in both cases (3 = 1, 4 the number 
Nf (2Nf) is even. According to the prescription of the case (3 = 4 we will do the replacement ./V — > 2JV 
at the end of the derivations. 

Following the recipe described in last section we now construct the kernels S^j\xi,Xj), j3 = 1,4, 



defined by the weight function (|39|). 

With the polynomials H^ f (x), given by ( pl| ) and ([42]), and ui^J(x), given by (|39|), we have 

<p i (x) = H^(x)J®(x)K (47) 

i = 0, 1, . . ., and Tt = span{c/?o, (pi, ... , ipN-i}. In each ensemble j3 = 1, 4 we now have the kernel 
K N (x,y), through ©. 

Because Nf is even in both ensembles we have J$.{x) = e - cx2+N s Xn \ x \. It follows then, that 

x — y xy x — y xy 

The polynomial i?^ 1 (x) does not have a constant term, see ((42|), and therefore the function ( |22"| ) 
vanishes. We have 



= aiV /°° dy fl^(y)^i(y)(^ + 2c) ^ exp(-q, 2 ) 

J —oo xy 

AT poo N 

= a N ^ dy H N N l (y) PN -2(y) y N * exp(-cy 2 ) = 0, (49) 

^ J —oo 

where pn-2(u) = H^ f _ l (y)/y. The function (^) equals 

(H^(y)) 2 (^- + 2c) j/"/ exp(-cy 2 ). (50) 
-oo xy 
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It follows from the orthonormality properties that the integral of the second term in (|50| ) equals 2ccin- 
The integral of the first term gives 

a N ^ d y ( H N f (y)f y^- 1 exp(- C y 2 ) 

X J — oo 

= « J v^^^ (< / (2/)) 2 exp(- Cy 2 )r oo - a^-L r dy y N f exp(-c y 2 )Y 

- /'°° % y N f[2H$(y)Hp(y) - 2yc(H% (y)) 2 ] exp(-cy 2 ) = 0, (51) 

£ J —oo 



X 

implying that 



i 



B N {x) = 2ca N = V2Nc2. (52) 
We then have the functions 



A(jr) = --V'(x) = --(2cx - = -ex + B N (x) = 2ca N = V2Nc^, 

2 2 x 2x 

C{x) = -^-Bjv-i = 2ca N = B N (x) = V2Nc^. (53) 
ajv-i 

Now we have the matrix in (p6[), making it possible to determine the matrix A after the choice of a 
number 2n of functions ipi(x), see (pop. 
The function 

N f 

— TW = — - 2 ^x, 54 

has a pole at x = and at x = oo, with orders uq = 1 and noo = 1 respectively. Thus we have 
n = n oo + n o = 2 meaning that we have to calculate two correction terms. 

We must find n = 2 linearly independent functions ipi,ip2 £ Ti an< i 2 linearly independent ifts , ^4 6 W -1- . 



These functions must be written as linear combinations of the functions (27) and (|28|), with only the 
case k = 0, and each fy, i = 1, 2, 3, 4, is therefore a linear combination of the four functions 

tpN-i(x), <Pn{x), , • (55) 

X X 

The function fN-ii^) is a polynomial of degree (N — 2) times uj^j(x) 1 ^ 2 and thus it follows immediately 
that 

Mx) = <PN-l(x), (56) 
In addition, with the help of orthogonality it is easily shown that 

Mx) = <Pn(*), ^^eW 1 (57) 

x 

The requirements are simply that ^3 ^4 ^ TL and V3, V^-'-V'o and ^3, ?/>4-Ly?7V-i (see (p9|)), which 
follows from orthogonality. 



With the choices (|5q ) and ([57]) we now determine the matrix A through (|26|), (|53| ) og (30). Inserting 
(H), (f§) and @, into p| gives 



(x,y) = -OAT C(p N (x)(p N -i(y)+Cip N -i(x)ip N (y) + — J -(tp N (x)ip N - 1 (y) + (p N -i(x)tp N (y)) : (58) 

2xy 
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A Comparison with ( |30|) gives the matrix elements of the 4 x 4-matrix A. We get 



A 



I 



\ 










-apfc 









N 



"f n 

"ON 9 U 



-aN~2 






\ 



/ 



/ 



V 



o 
o 

#a/2 



-(f) 










-(f) v W 



-(f) 1/2 c- 








N\l/2, 



-(f) 



which is symmetric and fulfills equation (31). 

Because of (S3) the 4 x 4-matrix B is given by the elements 



B 



12 



(eipufa), Bi 3 = {eipi,if) 3 ), B 



14 



(eipi,ip4) and 



#23 = (£^2,^3), #14 = (#1, V^), -B34 = (£^3,^4) 






(59) 



(60) 
(61) 



The generalized Hermite polynomials -ffj / (x) are either even or odd functions and the same is true for 
the square root of the weight function wj^(x) 1 / 2 (because Nj is even), and thus the functions ipi(x), 
i = 1,2,3,4, are even or odd. A general calculation of the inner product gives (sf,g) = 0, whenever 
/ and g have same parity. From this it follows that 



which together with 



B u = (eipi,ip4,) = B 23 = (eip 2 , ^3) 
results in 



/ 



B 





(#1,^2) 
(#1,^3) 







-(#2,^4) 






-(£^3,^4) 



0, 





(£^2,-04) 
(£^3,-04) 





(62) 



(63) 



We have 

Bij = (E1pi,1pj) 



dxipj(x) / dye(x-y)tpi(y) 



dy+ / dy)- 



from which it follows that 



dx iffj(x) / dy ipi{y), 



dy]tpi(y), 
(64) 

(65) 



whenever tpi is an even function. Notice that in every element, B^ ^ 0, one of the functions tpi or 
ipj is always even while the other is odd, which together with the feature B^ = —Bji implies that 
the absolute value of all matrix elements B^ can be written on the form (|65|). Leaving these matrix 

elements unsolved for a while, we now focus on the structure of the kernels S^\x, y), (3 = 1,4. 
Having the matrices A, B and (n = 2) 



J 



( 1 \ 

10 



V / 



C = J + BA, 



An 






N\l/2 



-(f) 



C2 









1 



-(f) 1/2 



(66) 



c 2 . 



and so on, the two kernels are given by 

4 



Sffix, y) = K N (x, y) - fo(x) Y,i A o C ooCohi eMv) ~ M%) Y\ A * C ™ C ^ £ Mv), (67) 



1=1 



k=l 
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S$(x,y) = K N (x,y) - - BAC)- 1 ^ e^v) ~ M*) " BACy l ] k2 eMv), 



1=1 



k=l 



where 



eipi(x) 



dy 



dyjipiiv) 



dy + dy 



dy)^i{y), 



1, 2, 3, 4. If ipi is an even function we see that 



while if ipi is odd then 



dyipi(y), 



dyipi{y). 



(68) 
(69) 

(70) 
(71) 



Now we turn to the troublesome integrals in the kernels. The explicit expressions for the four functions 
ipi(x) are given by <fj(x) = Hj 3 \x)u^}(x) 1 ^ 2 , (56) and (]57|) and are 



ip 3 (x) 



(-1)-- 1 

tp!(x) 



:f-l)! 



-IT 



^3 (a;) 



r ,7V , iV/+l 
L V 2 2 

r ( f + ^ti; 



1 3+N f ^£+1 

2 T 2 C 2\ — 

C * ^jvln ( C2; ) X 

2 1 



JV f 2 
2^+l e - C ^ 



c 4 L N 2 



(ex )x 2 e 



(72) 



For relatively large TV it is clear that an actual calculation of the two kernels sffi(x, y) becomes quite 
involved due to the presence of the integrals QS^), ( f70| ) and (71). Actually, it is the microscopic limit 
of the kernel 

(73) 



S { S\x,y) t 



l im -Ls W (— 



0=1,4, 



giving all the microscopic correlations, which has physical interest. Therefore it is necessary to deter- 
mine whether it is allowed to substitute an asymptotic relation for the polynomials in the integrals 
when the microscopic limit is taken. If the answer to this is positive, we can get rid of the iV de- 
pendence and derive a closed analytical expression for the microscopic kernel (|73|). With the help 
of Lebesgue's Majorant Theorem this question is pursued and answered in Appendix A. For (|70|), of 
course, the use of an asymptotic relation for ipi is allowed when the microscopic limit is taken, because 
the integral is finite for all N. We must check that the second equality in the following equation 



lim N X B. 



as well as the last equality in ( see uH])) 



lim N x 2 



dx ipj(x) / dy ipi{y) 



dx lim N x ipAx) / dy ipi(y), (74) 

N^oo 



lim N s e^(-) 

N^oo iv 



lim N d / dyiPi{y) 

N— >oo 



lim N d — 

N->oc N 



duipi 



N' 



du lim N S -^(-), 

N^oo iv 
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are valid. Here we assume that the power of iV A and N s have the exact values needed to make the 
number N^Bij and the function convergent for iV — » oo. This is a consequence of the 

assumption that the microscopic kernel (f73|) is finite. 

Now let us examine ([74]) for i = 1 and j = 2. After two substitutions we have 



foo JV /+ 1 , N f . rz J, / T± JV f 

B l2 = C N dz L N 2 _(2z)z2 { —-^e- z du L* (2u) U — e~ u , (76) 

JO 2 1 JO 2 1 



where all irrelevant factors have been collected in Cat. (We remember that (|6q ) only holds when ^ 
is an even function, meaning that Nf/2 is odd in case of i = 1 (see ( |72D ). But the fact .B^- = — -Bjj 
implies that all By can be represented on the form (|55|), when we keep track of the sign.) Using the 
result of Appendix A. 2 we find that ( |74| ) does not hold for B\2- From (122) and point (2) at the end 
of A. 2, this is a consequence of 

✓ JVf + 1 1 Nt\ 1 , . 

-(^ + 4"x) = -2^- 1 ' W 

and 

/N f + 1 1 .N f l x \ , . 

+ 7"(-/-o) )=-!>-!■ (78) 



V 4 4 v 4 T 
A similar treatment of the explicit expressions for the three other matrix elements -B13 , -B24 and 



gives that the corresponding equation (74) does not hold for any of them. An undesirable consequence 
is that we have to solve all four integrals B\2, -B13, -B24 and -B34 before enlarging N in the kernels (|67[) 
and ©. 

Turning to the question of the validity of the last equality in (|75| ) we must use the result of Appendix 
A.l. Consider for instance ipi and assume that this function is odd, that is ((Nf/2) + 1) is odd. Then 
from (J7D and © we have 



^ } = -iV/i dU ^N )cC -Nj x du Vi^FHl) 



e 



ex 



00 «/j 



2 " 



x (2t)t 4 e - { , (79) 



where all irrelevant factors have been skipped. From ( |109|) and ( |120| ) it now follows that ( |75| ) is noi 
valid for i = 1. This is caused by the fact 

N f + 1 ^Z.^ 3 (m 
" T ~ 4 " I" (80) 

Considering the functions eipi(x), i = 2,3,4, with ipi{x) odd (which of course is not possible for all 
functions at the same time), an analogue argumentation gives the same conclusion. When the function 
ipi(x), i = 1,2,3,4, is odd we therefore have to solve the integral eipi(x) to make an exact calculation 
of the kernels (p8|) and (pTj) possible (within reasonable time) for large N. 



We have learned that for finite N a calculation of -B12, -B13, -B24 and -B34 and an analytical expression 
for the four functions eipi(x), i = 1,2,3,4, are needed to make the results (|6S|) and (|67[) usable. This 
requirement has been met in Appendix B. 

In Appendix B.l we derive an expression for the function 

£ [aAn] (x)=e L«(x 2 ) X Pe-^, (81) 
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for integers n, /3 and a > — 1. The functions eij)i(x), i = 1,2,3,4, belong to this class of functions and 
it is easily derived that (see (|72|)) 



eipi(x) = k\c ±£rN f +i 



E1p 2 (x) : 

eip 3 (x) 



k 2 c±£. 



Nj_ v_ ll (c 2 s), 

2 > 2 ' 2 J 



fee 4 £r 



' 2 ' 2 



(C2X), 



with the coefficients 





= (-1)- 


k 2 


= fcl, 


h 


= (-1)- 


k-i 


= A*. 



(f-1)! 



(82) 

(83) 
(84) 
(85) 



r(f + 



:f)! 



r(f + 



(86) 



For each ^ the associated function <? above is given by ( |144| ), From (72) it is clear that if Nf/2 is 
even, then ipi and ip4 are odd functions and ip 2 and ip3 are even functions. For odd Nf/2, of course, 
the situation is the other way around. Thus for each kernel ( |68|) and (|67|) our derived expression will 
split into two parts, one for even Nf/2 and one for odd Nf/2. 
In Appendix B.2 the result for the function £ is used to calculate 



dx (j>Ax) / dy e(x - y)4>i{y) 



(87) 



where 



<f>j(x) = L^x^^-t og <j)i(x) = L%(x 2 )x p e 



3 - — 

V e 9. 



and n, m, (5, (3 are integers, fulfilling certain conditions. We have always an odd j3, meaning that (pi 
is odd. By keeping track of a sign we can always, with the help of B^ = —Bji as described earlier, 
secure that the function in the inner integral of the double integral Bij = (eipi,ipj) is odd (or even). 
Thus the four relevant numbers Bij are contained in the derived expressions for ( |87|) in Appendix B.2. 
Putting the functions (72) into (|87|), and remembering the mentioned possible change of sign when we 
secure that the odd function comes in the inner integral of Bij, gives the following result for the four 
matrix elements of B : 

Bn = (-l)-kjB 12 , (89) 



B\3 
B 2 4 '- 



(-1) 2 C 2klk 3 Bl 3 , 

N f , 

(-l)—+ 1 C 2k 2 k 4 B 24 , 



B^ = {-l)- + \k z fB, 



(90) 
(91) 
(92) 



where the numbers Bij are found with the help of table || and 
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even -f- 


given by eq. 


with parameters 


B12 


(|159|), without the (+)-term 


iVf+l a N f q Nf . -, W i 

a, a = — ^ — > P = 2' P = 2 + b n — m — — — 1 


B 13 


(fL59D, without the (+)-term 


a- V. «- JV/ 2 +i > )9 n-% 1 


24 






034 


(159), without the (+)-term 


a = a = — ^ — , p = — 1, p = -g-, m = n = y 



Table 2: A recipe giving Sy for euere Nf/2 . 



odd^ 


given by eq. 


with parameters 


012 


(|159j), without the (+)-term 


a, a = JV/ 2 +1 , /3 = ^ + 1, (3 = n — m = ^ — 1 


^13 





-59|), including the (+)-term 


a - iv { J +i , a - J %"\ /3 - 7 + 1, /3 - JV /, m - £ 1, n - £ 


024 


( 


159Q, without the (+)-term 


- O 2 7 — 2 ! — 2 ' — 2 ' ^ — 2 


034 


( 


159J), without the (+)-term 


J\f—l a N f 3 N, 1 AT 

a = a = — , P = -gS P = -j- — 1, m = n = — 



Table 3: A recipe giving 0y for odd Nf/2. 



Having calculated the four matrix elements of B and the four functions eipi we can now construct 
the two kernels (^38|) and (]67| ) of the non-xGaussian ensemble. The proceeding points summarize the 
construction of the kernels : 

(1) Through equation (|l"0|), the kernel K^{x,y) is given by (|47|) with the associated comments. 

(2) Assume that (-^-) is even (odd). Then it follows from (|72|), that ipi and ip^ are odd (even), and 
ip2 and V>3 are even (odd). The functions eipi, eip2, zips and eip4 are given by (|82|), (|83|), (|~~ 
and (|85[), respectively, with £ given by (P-44Q both for even and odd ipi. 



(3) The matrix elements -B12, -B13, -B24 and -B34 depend on whether Nf/2 is even or odd. These are 
given by @, @, (0) and ©, with B^ given by tables | and |. Matrix B is given by (§§). 



(4) The matrix 

C = J + BA 

is given by fl5S|), (|63|) and (|66|). We need the matrices 



AqC Q qCq for /3 = 4, and 
^C(I-SAC)" 1 for = 1. 



(93) 



(94) 



(5) Before putting all pieces together in the construction of (68) and ( |67| ) we remember the replace- 
ments 



(3 = 1 : 
(3 = 4 : 



c — > 2c and 
iV -> 2JV, 



2iV 



(95) 



in all parts above and in the functions (72). In the two ensembles (3 = 1, 4 we have c 



NT, 2 /3 



The five points above leads to two formulas for each sffl(x, y), (3 = 1,4 : one for even Nf and one for 
odd Nf for /? = 1, and one for even Nf/2 and one for odd Nf/2 for (3 = 4. 
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5 The microscopic spectral density 



In the microscopic limit the spectral correlations between the eigenvalues of the Dirac operator in 
the two Yang-Mills theories mentioned in the introduction can be derived from the corresponding 
symplectic or orthogonal non-xGaussian ensembles. In this section we present our results for our 
numerical expressions for the microscopic spectral density. 



From (14) it follows that the 1-point correlation function, that is the spectral density, is given by 

pW(x) = < } (x) = S<g\x,x), (3 = 1,4. (96) 

The machinery from the last section gives bothS { p(x,y), (3 = 1,4, and we therefore have in principle 
all correlation functions and especially the microscopic spectral density in the cases (3 = 1 and (3 = 4 
of the non-xGaussian ensemble. The microscopic spectral density is defined by 

^(,)Elim>(^ = 1,4, (97) 

N^oo is iv 

and this is of physical interest. We have not yet extracted the N dependence of our results for the 
Bij (derived in appendix B.2) and for the functions eipi (derived in appendix B.l) in the microscopic 
limit and we therefore have not an exact analytical result for ps(x), (3 = 1, 4. However, we can easily 
present plots of the scaled spectral density 

which, of course, converges towards pf\x) for large N. 

Due to the presence of the Vandermonde determinant in the probability distribution one expects that 
the spectrum becomes more rigid when the parameter (3 is increased (see (§)). Thus we naturally 
expect the scaled spectral density p$(x) in the ensemble non-xGOE to be flat compared to p$ in 
non-xGUE, and that p$ in non-xGSE is the most oscillating function of the three of them. These 
features are of course valid for all values of N. 

In both ensembles (3 = 1, 4 the kernel sffi (x, y) is equal to a kernel Kr 2 )N ( x > V) from the unitary ensem- 
ble (3 = 2 (remembering the modifications associated to each value of (3= 1,4), plus n corrections, in 
general (see ( |T5| ) and (|16D). Thus in the two non-x ensembles we have p^\x) = p^\x) + P^corv 
(3 = 1,4, where p^\x) = (N~ 1 )K( 2 )n(N~ 1 x, N~ 1 x) and Pn\ ovi (x) constitute the 2 correction 
terms. When having calculated the matrix B and the four functions eipi{x) both scaled corrections 
/°7V^corr ( a; )' P = 1> 4, are given by the two last terms on the right hand side of ( |67| ) and (^). The con- 

(2) 

tribution from the scaled spectral density p N (x) is an oscillating function in both ensembles (3 = 1,4, 

and the flat spectrum of p$ (x) in the ensemble non-xGOE is therefore expected to be a result of the 

fact that py- corr (x) precisely cuts off the peaks of p^\x), and at the same time the term p^ corr (x) 

must cause a highly oscillating pff(x). When plotting p^\x) together with the associated p^\x) we 
indeed observe these highly non-trivial features. The results are illustrated in figure |l| and |. 

In the limit N — > oo we have pffl(x) — ► pi^\x), (3 = 1,4. This convergence, however, is especially fast 

for x close to zero. This is illustrated in figure^, where we have plotted p^\x), (3 = 1,4, for different 
values of N. The spectral density is an even function (see (|40|)) and for finite N we have the so-called 
half circles (with a hole in the vicinity of x = 0, due to the term x Nf in the weight function), which 
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Figure 1: In non-xGOE we have p N (x) = p N (x) + p N corr .( x )i where p N (x) is from the associated 
unitary ensemble (given by orthonormal polynomials with respect to uj^ ; {x) 2 ) and p$ ' corr (x) is the 

(2) 

scaled correction derived in last section. In left column the oscillating term p\ (x) and the entire 
P4o( x ) are plotted together for N = 40 and the two values Nf = 2 (above) and Nf = 4 (under). 
The corresponding correction terms P40 cott.( x ) are plotted in right column. We observe that the flat 
spectrum of non-xGOE is exactly a result of p^ corr _(x) , cutting of the peaks of p^(x). In addition 
the curve of p^J (x) moves away from zero when Nf is enlarged, due to the presence of the determinant 
function in the probability distribution (see (^)). 
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. 4 



0.3 



0.2 



0.1 





Figure 2: In non-^GSE we have p$(x) 



corr where p\y(x) is from the associated 

,( 4 )/ 



unitary ensemble (given by orthonormal polynomials with respect to oj n {x) and the replacement 

N — > 2N) and pj\^corr.( x ) ^ s the scaled correction derived in last section. On the left the term p^{x) 
and the entire p^(x) are plotted together for N = 40 and the value Nj = 4. The corresponding 
correction term P40 C orr. ( x ) * s plotted on the right. We observe that P40 corr.( x ) gi ves a highly 
oscillating p^(x) (a rigid spectrum). 



is clearly seen in figure [3|. Because of the the Gaussian distribution the microscopic spectral density 
fulfills p s (x) — > 7T _1 for x —>■ 00 and thus the curves p^ (x), (3 = 1,4, are of order ~ 7r _1 before the 
half circles start approaching zero (see also figure [l] and ^) . 
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Figure 3: The scaled spectral density p N (x), for increasing N, on the left for /3 = 1 and on the right 
for (3 = 4. The convergens in the vicinity of x = is clearly observed and we have p^ \x) ~ pi (x) 
near z = as a consequence. In the limit iV ^ 00 we have pff (x) — ► ps (x) and in the limit x —> 00 
we have pf\x) — ► ir^ 1 . 
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5.1 The first spectral sum rule 



The spectral sum rules derived in field theory are all given by the microscopic spectral density p s of 
the Dirac operator [27j and the result for p s derived from RMT can therefore be verified through 
the knowledge of a spectral sum rule. The first sum involves the sum over the eigenvalues of the Dirac 
operator in second inverse power and in the limit N — > oo we have trivially 



( 



1 



00 dz p M. 



(99) 



where the sum is over all Xk 7^ (both positive and negative) and the average is taken over all gauge 
fields. For both field theories corresponding to the ensembles = 1,4 in (2 + 1) dimensions (see table 
|l|) the same spectral sum rule has been derived in [10, [ll]] an d it reads 



(E 



1 



4 

w 



(100) 



According to this, the first spectral sum rule involving the eigenvalues of the Dirac operator in two 
completely different field theories in (2 + 1) dimensions should be identical. 



With the help of our numerical expressions for p, 



(0), 



for the microscopic spectral density p, 



O 9 ), 



x) m 



the two ensembles = 1,4 we thus have a consistency check on our results through (|99| ) and (|100[) . 
From figure ^ and |3] it is easily seen that the effect of the finite iV is very small for a relatively large 
N. As seen our results does not agree with the spectral sum rule ( 100 ). In the ensemble = 4 it 



N f 


f 00 J-v PiO ( z > 
J-00 az z 1 


(4), v 


4 

N f 


2 


0.402 


0.946 


2 


4 


0.148 


0.436 


1 


6 


0.092 


0.297 


0.667 


8 


0.067 


0.227 


0.500 


12 


0.044 


0.155 


0.333 



Table 4: The calculated first spectral sum rule, i.e the right hand side of 
= 1 and = 4 together with the spectral sum rule (llOi 



in the two ensembles 



seems the spectral sum rule approaches 2/Nf for large Nf. 

Next we turn to the most elementary check of our results, which in fact is an actual generation of 
matrices. 

6 The Monte Carlo simulation 

We have performed a Monte Carlo simulation in order to numerically verify the derived scaled spectral 
densities p^\x), = 1,4 The matrices in the generated ensemble are distributed according to (|3|), 
with all m,f = 0. The simulations are done for various values of the matrix size, N, and different 
values of Nf in the two cases = 1,4. 

We diagonalize matrices from an ensemble made by a simple Metropolis algorithm. The eigenvalues, 
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Figure 4: Scaled spectral density of the non-xGOE ensemble, compared with the Monte Carlo data, 
N = 12, Nf = 2 and 4. 

x, are then collected to form a histogram, h(x), which is rescaled 

x — > xiV = z, 

h - | (101) 

to obtain jjh(4j). The rescaled histogram is then compared to the plot of the analytically obtained 
scaled spectral density. Note that we compare the analytic plots for finite N with the corresponding 
simulation data. We thus expect the plots to match for the whole spectrum. For large N the scaled 
spectral density close to z = has the form of the microscopic spectral density. 
We handle the two cases (3 = 1,4 individually as follows. 

6.1 13 = 1 

The Metropolis step consist in changing the entries in the matrix and accepting these changes according 
to the action 

S(M) = — TrM 2 - N f In ( det M) . (102) 

The changes is done in a way that explicitly maintain the symmetry of the matrix (M = M T ). 
In this way we obtained an ensemble of 100,000 effectively uncorrelated matrices for N = 12 and 
10,000 matrices for N = 40, both for Nf = 2 and 4. The resulting histograms together with their 
analytic equivalents can be seen in figure |4| and ||. We observe excellent agreement in all parts of the 
spectrum, for all values of the parameters N and Nf. 

6.2 = 4 

In this case the Metropolis steps are accepted/rejected according to 

S(M) = NTiM 2 -Nfhx{ det M) (103) 
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Figure 5: Scaled spectral density of the non-xGOE ensemble, compared with the Monte Carlo data, 
N = 40, N f = 2 and 4. 

A matrix from the symplectic ensemble can be represented by a 2N x 2N matrix, M. This matrix 
can be written as a sum of direct-product matrices 

M = M e + Miei + M 2 e 2 + M 3 e 3 (104) 

where 

- = (;!)■ <**> 

and the coefficients Mi , M 2 , M3 are real antisymmetric N x N matrices and Mq is a real symmetric 
N x N matrix. In this way the matrix M has 2N eigenvalues which are doubly degenerate. 
The Metropolis update works on the four real matrices, Mi, and M is then constructed from these 
and diagonalized. In this the way the symmetry and structure of the matrices is kept at all times. For 
the histogram only the iV different eigenvalues is used. 

When working in this 2N x 2N representation the action, ( |103| ), has to be changed. The trace gets 
twice as large and the determinant is squared in this representation, so an overall factor of 1/2 should 
be included: 

S(M) = — TrM 2 - -J- In ( det M) (107) 

These simulations are about twice as slow, so the maximum matrix size is half as big. We have done 
runs for N = 8 and 20, with Nf = 2 and 4, with the same statistics as in the (3=1 case. The plots 
are shown in figure || and |^. Here the agreement is even more striking. Zooming in on the first few 
eigenvalues in the iV = 20 case we see absolutely perfect agreement indeed, see figure || 
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Figure 6: Scaled spectral density of the non-xGSE ensemble, compared with the Monte Carlo data, 
N = 8, N f = 2,4. 
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Figure 7: Scaled spectral density of the non-^GSE ensemble, compared with the Monte Carlo data, 
N = 20, N f = 2,4. 
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Figure 8: Same plot as figure [?] but zoomed in on small values of 



7 Conclusions 



The purpose of this study has been to derive the kernels (x,y), (3 = 1,4, from which all spec- 
tral correlation functions of the non-xGaussian orthogonal (non-xGOE) and symplectic (non-xGSE) 
ensembles can be determined. In the microscopic limit the two ensembles describes the microscopic 
spectral correlations of the low energy Dirac operator in two Yang-Mills theory in (2 + 1) dimensions. 
Thus the universality class non-xGOE carries information about the spectral correlations in a SU(2) 
gauge theory with fundamental fermions, while the universality class non-xGSE does the same for the 
SU (N c ) gauge theory, with N c arbitrary, and fermions in the adjoint representation. 
With the help of Widom's new method we have derived the kernels <sjy (x,y), in non-xGOE ((3 = 1) 
and in non-xGSE {[3 = 4) for massless fermions. For non-xGOE our result are valid for all inte- 
gers Nf and for non-xGSE we have a result for even Nf. We plotted the scaled spectral density 



Pn \ x ) = N l x,N for different values of the parameters N and Nf. The plots of 

Px{x) possesses the expected features, a flat spectrum for j3 = 1 and a highly oscillation spectrum 
for (3 = 4, and indeed we see perfect agreement with our computer simulated spectra. For large N the 
scaled spectral density p^\x) coincides with the microscopic spectral density p s {x) in the vicinity of 
x = 0. Our calculated spectral sum rules from p^\x) of course match the computer generated ones 



perfectly. We have not found agreement with the sum rules of [[H], 11]. 



Recently the microscopic spectral densities of the massive non-xGOE and non-xGSE have been de- 
rived by a different method p3] . In the massless limit these result seem to match ours for large N. 

Acknowledgements : Poul Henrik Damgaard is thanked for many fruitful discussions. The work of 
C. Hilmoine is partially supported by a Rosenfeld Fellowship. 
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A Interchange of integration and the microscopic limit 



In this Appendix we look at the large n-behavior of two types of general integrals involving generalized 
Laguerre polynomials. Specifically we examine when interchange of integration and the microscopic 
limit is allowed. To this end Lebesgue's Majorant Theorem is used. The first section is relevant for the 
calculation of erpi in the limit N — > oo, see (|75|). The result of section 2 is relevant for the calculation 



of the elements of the matrix B in the limit N —* oo, see (74). 
A.l Jim / dy n il+a) - a L%{2y)y a exp(-y) 

•J 

We define 

poo 

T% a =\ dyn^~ a L^2y)y a exp(-y), (108) 
J o 

where n is a positive integer, a > — 1 and a are real numbers. We wish to calculate 

limJ°' a =lim / dy n {1+a ^- a L°;(2y)y a eM-y), (109) 

n— >oo n-^oo Jq 

and therefore examine when the interchange of the limit and the integration is allowed. In order to 
do this we employ Lebesgue's Majorant Theorem. At first we convince ourselfs that the sequence 
of functions in the integrand is convergent. The integral in ( |109| ) is rewritten by the substitution 
u = (2ny) 1 / 2 , and we get 

/•oo „,2 „.1 

lim X£' a = lim / dun ( - 1+a ^- a L^( — )u 2a+1 n- a - 1 2- a exp(- — ). (110) 
n^oo n n-^ooJ Q n n 2n 

The integrand in ( [110D 

f n {u) = n^- a L a n {-)u 2a+l n- a - l 2' a ex V (-^-l (111) 
n 2n 



is convergent for n — » oo : Inserting the well-known asymptotic formula [28] 



lim -L L %(-) = x-<*J a (2x), (112) 
n— >oo n a n 



2 2 

lim Uu) = lim n( 1+a )- a L"( — )u 2a+1 n- a - 1 2- a exp(- — ) = 2- a u 2a+1 - a J a (2u). (113) 

n— >oo n^oo 71 2n 



we get 



Thus, the integrand in ( |109[ ) is convergent. We must now look for an integrable majorant M(y), 
fulfilling for all n 

\Uv)\ = \n^- a L«(2y)y a eM-y)\ < M(y). (114) 
If we can find such a function A4(y), the Lebesgue Theorem states that interchange of the limit and 



integration is allowed in ( 109 ). In this case we can insert (113) in ( |109| ). 

We start by splitting up the integral into integrals, in the first integrating over the interval ]0, e[ for 
some finite e and in the second over ]e, oof. Since the first integral is finite for every n we focus only 
on the second integral. By choosing a sufficiently large n' we can use an asymptotic expression for 
the integrand. We then check if the absolute value of this expression has an integrable majorant for 
all n > n', see (114). Using the asymptotic property (valid for x > 0) [pq] 

rv/ \ 1 2. ( 2a-l \ ( 2a-\ \ 1 07T 7T . 2q-3 , 

L n\ x ) = — e2x 4 n 4 J cos[2(nx)2 —] + 0(n * ). (115) 

7T 2 4 
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for L%(2y) in (£6g), we get 



lim T% a 



lira 

n— »oo 



dy n (1+a " a) 

x |vr" 1/2 exp(y)y" (2£ ? i) 2" ( 
n -l/2 2 -^) Hm n (l+a-f-|) 



2a + 1 
4 



'n v 4 ' cos 



2(2ny) 



a7r 
~2~ 



}y a exp(-y) 



2a + l 

y 4 " 



cos 



2(2ny)5 - — - - 



(116) 



for n > n' . 

We know that the sequence of functions in (jH6|) is convergent, and the function 



M(y) 



C IT 



-l/2 9 -( 



2"l — ) y" 



c> 1, 



(117) 



is then a majorant to the absolute value of the integrand in ( |109| ) for n> n' and some c > 1, which take 
account of the term O in (|115| ). Since we need / e °° dy M.{y) < oo we must have — (a/2 + 1/4) + a < — 1. 
So we conclude, if 

a/2 -a > 3/4, (118) 



i/ien ( |117| ) is a integrable majorant and we can then interchange the order of the limit and the integral 
in ( \10Q ). Asumming this in (113), gives us 



lim T°' a 

n—>oc " 



du U 



2a+l-a 



(119) 



for a/2 - a > 3/4. 

On the contrary, if a/2 — a < 3/4, then the integrand in ( 109 ) has no majorant, «M(y), which fulfills 
Jq 00 dyM{y) < oo. This follows from the fact that the integral from e to infinite of the absolut value 
of the integrand in (|116| ) is not finite. Substituting z = 2ny in ( 116 ) gives an integrand on the 
form z~ s cos(2z 1//2 + fe), where 5 = —(a/2 + 1/4) + a (the n dependence vanishes of course). The 
smallest majorant to \z~ s cos(2z 1 / 2 + k)\, is C\z~ s cos(2z 1 ^ 2 + fc)|, for some C > 1. This majorant is 
not integrable when 5 = —(a/2 + 1/4) + a > —1, and thus we have no integrable majorant. When 
the absolut value of the asymptotic function does not integrate to a finite number, then of course no 
majorant exist for which the integral is finite (see ( |114j) ). 
We conclude, when 

a/2 - a < 3/4, (120) 



interchange of integration and the limit in ( |109D is not legal. In this case we first have to solve the 
integral in ( |109j ) and subsequently derive an expression in the limit n — > oo. 



A.2 



lim 

n— »oo 



dt n A L^(2t)^exp(-t) / du L*(2u)u p exp(-u 



We will now look at the limit n — > oo of integrals of the following type 



T 



■a,p,a,P 



n 



dt L^(2t)t /3 exp(-t) / du Z£(2u)i/exp(-u), 



(121) 



where a, a,/3,/3 > —1 and A is real and n is a positive integer. The matrix B in section 5 is given by 
integrals of this type, and we are interested in an expression for B in the limit n — > oo. It is assumed 
that the integrand is convergent for n — > oo, meaning that A has a certain value, the size of which is 
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irrelevant for the question addressed in this Appendix. 

Analogous to the previous section we investigate whether the limit and the integration can be inter- 
changed in the expression 



lim T, 

n—too 



roo r-t - 

lim / dt n A L^(2t)t /3 exp(-t) / du L^(2u)u^ exp(-u) 
n ^°°Jo Jo 



By assumtion the sequence 



(122) 



/„(*) = n x L%(2t)tP exp(-t) / du L*(2u)vP exp(-u), 

Jo 



(123) 



is convergent ( compare with ( |111| ) og (113) ). We must determine whether or not the absolut value 
of the integrandt in ( [L22|) , that is \f n \, has an integrabel majorant «M(t). Like in Appendix B.l we 
persuit this question by choosing a large n and put in the asymptotic relation ( [11 5| ) for both Laguerre 
polynomials : 



lim 2?'^ » 



lim n 



A 



dt \ U2 4 t ^2+4 Z 3 ) 



cos 



■ . i air 7T' 
2(2nt)* - — - - 



du 4u (2 + 4 z 3 ) 



cos 



' . 1 «7r 7T 

2(2nu)2 

. v 7 2 4. 



(124) 



Like in Appendix B.l we skip all contributions involving irrelevant integrals over ]0, e[, because in 
these cases the exact integrals are finite. The sequence in ( |122| ), that is f n above, is convergent by 



assumtion. Thus for all n the absolut value of the integrand in (124) is always smaller than the 
function 

(125) 



where c > 1. Letting 5 



1 "— '" 1 1 (3) we have trivially 



-(f + andp 



'a 1 1 
- 2 4 



dt r 



du 



We see immediately that for 
P + 5 + 1 



a + a — 1 



+ /3 + /3<-l and 5 



•(f + i-ffl<-i 



(126) 



(127) 



the integral in ( |126| ) is finite. For large n we therefore have that the absolute value of the integrand in 
( p_23|) always is smaller than the function ( |125| ), and if (127) is valid the function ( |125| ) is an integrable 
majorant. 



On the contrary we see from ( |124| ) in case of 5 = — (^ + \ — /3) > — 1 and p = — (^ + | — p) > —1, 
then the integral over t is divergent. It follows that the absolut value of the integrand in ( |122[) does 
not have an integrable majorant, since the asymptotic expression does not even have one. 
The other cases of values of p and 5 are not relevant for us and we ignore them. 



Summarizing the discussion we have 



(1) I f (p.27| ) is valid, then the absolut value of the integrand in (122) has an integrable majorant 
( |125| ). Thus interchange of the limit and integration in (122) is legal: 



lim Z?,/W 



00 rt 

dt lim n A L°(2t)t /3 exp(-t) / du L™(2u)u p exp(-u) 
n -*°° " Jo 



(128) 
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(2) If -(§ + \ - p) > -1 and -(§ + \ - p) > -1, then the absolut value of the integrand in Q 
has no majorant -M(i) fulfilling f£°dt M(t) < oo. Therefore interchange of the limit and the 
integration in (|122[) is not allowed. In this case we must solve the integral before the limit is 
taken. 



B Two integrals solved 

In this Appendix we find expressions for two different types of integrals. We wish to derive an 
expression for the function 

™2 poo _ - ~,2 

£[a,An](*) = * ^Vexp(-y) = j Jy £ {x- y)L«(y 2 )/exp(-|-), (129) 

where a > — 1 is a real number and n, (3 are integers. 

For integers m and n we calculate the following numbers 

™2 2 roo 2 -2 

Bij = (e L^x 2 )^ exp(~) , (aVe"^) = / dx L^{x 2 )x p e L^x 2 )^"^ 

& J — oo 



oo 

OO 2 /"OO _ ,,2 



OO <7 — oo 



dx Z4(xVe"^ / dy e(z - l/)i£(j/Ve~ V , (130) 



The index i refers to m, a, and /?, while j refers to n, a and /?. m, n, (3 and /3 are integers and a and 
a are real numbers greater than -1. 

The function (|12S| ) is a part of the inner product in Bij and first we therefore derive an expression for 
L29D. 



B.l £ L«(xVexp(-y) 
We start by defining 

E(\,x) = ex x (xe~^y (131) 
where A > —1 is an integer. We have the following recursion formula 



E(X,x) = / dy e{x-y)y x (ye ^) 

J —oo ^ ' 

■ f x f°° 1 \( 

/ dy - / dy y [ye i ) 



12 12 /*00 , 2 

,.A- 2 7 



2 /■OO 

-e 2"x A + A / dy e(x — y) y 

J —oo 



ye 



2 



= -e~^x x + XE(X-2,x), (132) 
for A > 1. For A = —1 we have 



roo 

E(—l,x) = / dy e(x — y)e 2 

j —00 

- ^{/:>/>-f } 



e 2 



do 

X 2 







r/,/r-~ - \/Jerf(-|). (133) 
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-e 2 



And for A = we get 

roo _y2 

E(0,x) = / dy e(x — y)ye 2 
J — 00 

For A even we immediately get from ( |132| ) and ( |134D 

E(X, x) = -e~ir f x x + \x x ~ 2 + . . . + A(A - 2) • • • 4 • 2) . 



In the case of odd A, ( 132 ) and (|133j ) give us : 
£(A,x) 



-e 2 



(x A + Ax A " 2 + . . . + A(A - 2) • • • 5 • 3xJ + A(A - 2) • • • 5 • 3^|erf(-^=). 



Writing out the terms of the Laguerre polynomials through [25 

/ n + a\x m A , ^ m T{n + a + l) 



W = E ("!)' 



m=0 



n — m 1 ml 



E(-!) f 



m=0 



r(n — m + l)r(m + a + 1) m\ 



(134) 



(135) 



(136) 



(137) 



we get 



e L^(a; 2 )x /3 exp(- 



dy e{x -y)L^(y 2 )y f3 e 



n f°° r ( v 2 \ 

J2 a i d v £ ( x - v) y (ye'^ ) 

■ n J —OO V ' 



i=0 



a 



dy e(x - y) / 1 (ye ^ ) + a x \ dy e(x - y) y^ +1 (ye "2 \ + 
^ ' J — 00 ^ ' 

dye(x-y)y 2n+ P- 1 (ye- V 4), (138) 



. . . + a. 



where ai are the coefficients to x l in ( 137 ), 

Assuming that (3 is odd, all powers in (138) are even and we use ( |135| ) on each term to give us 



e -L"(x 2 )ar exp(- 



x 



-e 2 



a (x^ 1 + (p- l)xP~ 3 + . . . + 03 - l)(p - 3) • • ■ 4 • 2) + 

01 (x~ p+l + {j3 + l)^" 1 + + - 1) • • • 4 • 2) + 

+a n (x^ +2n ~ l + {fl + 2n - l^" 3 + . . . + (fl + 2n - 1)03 + 2n - 3) • • • 4 • 2 



-e 2 



{a + (/? + l)oi + 09 + 3)03 + l)a 2 + • • • + {fi + 2n - l)(/3 + 2n - 3) • • • (fi + l)a n } x 
(x?- 1 + (/3 - l)xP~ 3 + (p-l)(j3- 3)x^" 5 + • • • + (p - l)(/3 - 3) • • • 4 • 2) 
+^ +1 (a! + (/? + 3)a 2 + (^ + 5)(^ + 3)a 3 + . . . + {fi + In - + 2n - 3) • • • [fi + 3)a n ) + 
P +3 fa 2 + + 5)a 3 + . . . + (0 + 2n - 1) • • • {fl + 5)a n ) + ... 



... + X P +2n - l a r , 



(139) 
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When (3 is even we use (|13q ) in ( J138|) , and we get 



e L^(x 2 )x /3 exp(-^-) 



-e 2 



ao^" 1 + (/3 - l)x^" 3 + . . . + (/? - 1)03 - 3) • • ■ 5 • 3 erf(^=)) + 
ax + (0 + l) X P- 1 + . . . + (/3 + 1)09 - 1) • • • 5 • 3 erf(^) + • • 



+a n (x~ p+2n ~ l + ((3 + 2n- l)x^~ 3 + ... + (/? + 2n - l)(/3 + 2n - 3) • • • 5 • 3 erf( 



Using the definitions 



(140) 



(2k)\\ = 2k(2k - 2) ■ ■ ■ 4 • 2 = 2 fe T(A; + 1), 



(141) 



(2k - 1)!! = (2k - l){2k - 3) • • • 3 • 1 = vr"^ 2 fc T(k + - 



(142) 



where k is a positive integer, we can write (|139| ) and (|140| ) in a more compact way. 

For every a^, < i < n, the righthand side of eq. (|13S| ) is equal to — e~ x I 2 times a polynomial of 

order (/3 + 2i — 1). This polynomial can be rewritten as 



r(^ + ») 



(143) 



A similar expression exist for the polynomial in the framed part of ( tl40[) . Using these expressions for 
the polynomial parts of (140) and (139) and inserting the Laguerre coefficients of x l in ( |137| ) in 
( |140|) and (jl3|) leads to 



£ [a,p,n]W = £ L n(^ 2 )^exp(-y) = 
(-If r(n + l + a) 



-e 2 



[E 



i+i {-§} 



, :Q T(i + 1) T(n + 1 - z)r(a + 1 + 1) 
™ (_!)< T(n + l + a) 



x (/3-l-2j+2i) 2 j 



j=0 



r(*£ + i-i) 



^ r(* + 1) T(n + 1 - i)r(a + 1 + *) 



2 («-!)r(< + ^±i 



e 2 



(144) 



For even j3 we must include the terms in the brackets {. . .}, while these are neglected for odd (5. 
For (3 = a modification of the solution with the brackets is valid : for i = we must delete the term 
7r. All contributions from a negative upper limit in the summation are set to zero (for 



-xe 



-x 2 /2 



instance for (3 = 2 and i = the term involving 2~2j=o equals zero). 



30 



/OO 2-2 
dx L^(x 2 )x p e-^e L«{x 2 )x p e~^ 
-oo 

Using ( 144 ) we are now able to calculate ( p.30[) . Since j3, (3 are integers, the interior and exterior 
functions in ( |13C| ) 

f int (x) = L«(x 2 )xPe-£ 

and 

2 

f ext {x) = V^{x 2 )x^&~^ 

are either even or odd. In this case it is easily shown that Bij is non-zero only when (3 is even and 
(3 is odd or vice visa. In the following we assumed that (3 and (3 fulfills this. Using the relation 
Bij = —Bji we can make sure that the function in the interior integral in ( |130| ) is odd. If (3 is odd this 
is automatically the case and we have 



B^ = I dx fext(x)ef int (x), (145) 

J — oo 

with efi n t{y) given by ( |144| ) . In case of (3 even we simply use Bij = —Bji to switch the two functions 

/oo 
dx fint(x)ef ext {x). (146) 
-oo 

where now ef ex t(x) is given by (144). 

So from now on we assume (3 is odd and (3 even. 

Denoting by p ( g_ 1+2n the {(3 — 1 + 2n) order polynomial from ( |144| ) , we want to calculate 



/oo 
dx L a m {x 2 )x?e- X Pp_ l+2n (x) (147) 
-oo 

X 2 

Since both pp_ 1+2 n an d L%[x 2 )x^e~~ are even, we have 

roo „ 

Bij = -2 dx L^{x 2 )x p e- X p ^ 1+2n (x 2 ). (148) 

JO 2 

Making the substitution z = x 2 , gives 

f°° 0-1 _ i 

B^ = - dz L^(z)z 2 e z P0„ 1+2n (z2). (149) 

J 



Let us rewrite Pp_i +2n using (|139| ) 

1 0+2n-l 0+2ra-3 + 3 0+1 

Pp-l+2n( z2 ) = h p+2n-l z 2 + b 0+2n-3) z 2 + ■ ■ ■ + b (p+3) z 2 + b ((3+l) z 2 + 

0-1 0-3 . 0-5 



T(z— + b ( ^ 3) z— + b { p_ 5)Z — + • • • + bo) (150) 



The coefficients bi are given by products of (3 and dj, and T is the factor in brackets ({•••}) in ( 139 ). 
Using (|141|) and (|142|) T can be reduced to : 



{a + 09 + + G9 + 3)(0 + l)a 2 + . . . + (/3 + 2n - + 2n - 3) • • • (/3 + l)a n } 
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b (f3+2n-l) zn 1 + b(0 + 2n-3) z71 2 + • • • + fo (/3+3) Z + 



By insertion of (|150| ) in (|149| ) we get 

Z" 00 _ 3+3 

iBjj = - / ffc L*(z)e z z 2 
J o 

-T ^ L* (z)e-^— + b w _ 3) z— + 6 (/ g_ 5) ^— + . • • + bo) . (152) 

We restrict the parameters to the following four cases 

a = ~ ^ ~ > an d (n — 1) < m or (n — 1) = m, (153) 

/3 + /3 

a = — h 1, and (n — 1) < m or (n — 1) = m. (154) 

In case (|153j ), orthonormality reduces (|152|) to 

roo 3—1 / 3—1 3— 3 3— 5 \ 

Bij = -T J L^(z)e- z z— (z— + 6 (/ 3_ 3) z— + 6 (/ 3_ 5) z— + . . . + b ) , (155) 
for (n — 1) < m, while for (n — 1) = m we have 

/•oo 

Bij = —b(p + 2n-i) / ^ L m (z)e z z z — 

J 

f 00 / /3-1 5-3 3-5 \ 

T dz L^{z)e- Z z a (z— +b {0 _ z) z— +b Cp _ 5) z— + ... + b a ). (156) 
The first term in (|155| ) reduces to 

/■oo /-oo . ™| N 

-6^ +2 n-i)y o «te L-(z)e- 2 ^z m = -6 {/ 3 +2n _ 1} ^ cte L« ( Z ) e -* 2 «(^_L« (*)) 

m! m! r(a + l + ?n) 

= -ft (^+2n-l) /-Qro' 1 ™ = _6 (0+2n-l)(_-g m ~\ 

= -6 (/ 3 +2n _ 1) (-l) m r(a + 1 + m) = -a n (-l) m T(a + 1 + m) 

= -^(-l) m r(a + 1 + m) = - ti^(-i)»T(a + 1 + m) = ^ ± ^r^ - (157) 
n! (m + lj! 1 (m + 2) 



Working out explicitly what the 6j's are in terms of T-functions and using [28] 



dx e~*x^L»(x) = r[ r 7]r[1 1 + r /i + n " 7 j . (158) 
o nV ' Y[n + l]r[l + n - 7 ] V ; 



on each term in ( |155| ) gives us the result 
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n 



B 



T{S+l + i) T(n + l + a) 



i _V 2 

r(M)r(i + i)r(n + i-i)r(a + i + 

/3-1 



1 ' T(a-j)Y{l + m + j) , T(^±i) .r(a + m + l 



E 



r(m + i)^ r(i + j) 



2' 



+ 



r^-^" r(m + 2) ' 



for 



a, (n — 1) < m, and /? odd, 



(159) 
(160) 



The last term is only added in the case (n — 1) = m. 



In the case ( |154f ), only one term in the first line of ( |152| ) survives because of orthonormality. It 
reads 



6 (/ 3+i) / dz L a m (z)e z z ? = -b (B+1) \ dz L^(z)e z z 2 + z 



(/Hi) 



poo 

-b {m / dz L^ n (z)e- Z z a z- 1 = -b (m T(a). (161) 
J 



By (p9D and (]150D we have 



b {0+1) = (ai + {13 + 3)a 2 + {p + 5)(/3 + 3)a 3 + . . . + (/3 + 2n - l)(/3 + In - 3) • • • (p + 3)a n 
T - a 



(162) 



The contribution from the second line in (152) is given by (159) without the (+) term and with the 
substitution a — ► a — 1. Collecting the parts we get that for 















2 


— + 1 = q, (n — 1) < m , and /3 odd 




(163) 


we have 










By = 


Righthand side of eq. 


(159) with a —* a — 1 and without the term (+) 


+ 1 


r(a), 








(164) 


where T is given by eq. (|151|). 
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